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Abstract. The present status of the longitudinal polarized parton densities (PDFs) and the contri- 
bution of their first moments to the nucleon spin is discussed. Special attention is paid to the role of 
higher twist effects in determining the PDFs and to the polarized strange quark and gluon densities, 
which are still not well determined from the present data. 
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INTRODUCTION AND METHODS OF ANALYSIS 



Our present knowledge about the nucleon spin structure comes mainly from polarized 
inclusive and semi-inclusive DIS experiments at SLAC, CERN, DESY and JLab, po- 
larized proton-proton collisions at RHIC and polarized photoproduction experiments. 
One of the important and best studied aspects of this knowledge is the determination 
of the longitudinal polarized parton densities in QCD and their first moments, which 
correspond to the spins carried by the quarks and gluons in the nucleon. 

There are some peculiarities which make the QCD analysis of the data more compli- 
cated and difficult than that in the unpolarized case: 

i) A half of the present polarized DIS data are at moderate Q 2 and W 2 (Q 2 ~ 
1—4 GeV 2 , 4 GeV 2 < W 2 < 10 GeV 2 ), or in the so-called preasymptotic region. So, 
in contrast to the unpolarized case, the 1 /Q 2 terms (kinematic - y 2 factor, target mass 
corrections, and dynamic - higher twist corrections to the spin structure function g\) have 
to be accounted for in the analysis in order to determine correctly the polarized PDFs. 

ii) Due to the lack of charged current neutrino data, only the sum (Aq + Aq) can 
be determined from inclusive DIS and information from semi-inclusive DIS (SIDIS) is 
needed for the flavor decomposition of the polarized sea (Aw, Ad, As). Such an analysis, 
however, requires knowledge of the fragmentation functions which are not well known 
at present. Note that in the unpolarized case, SIDIS processes have never been used for 
a flavor separation of the parton densities. 

The best manner to determine the polarized PDFs is to perform a QCD fit to the 
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data on g\/F\, which can be obtained if both An and A± asymmetries are measured. 

In the case if only An is measured, the quantity An /D{1 + y 2 ) is a good approximation 

of g\/F\, where D is the depolarization factor and y 2 = 4M 2 x 2 /Q 2 . The data on the 
photon-nucleon asymmetry A\ are not suitable for the determination of PDFs because 
the structure function gi is not well known in QCD and the approximation [A\) theor [= 
g\/F\ — y 2 g2 /F\)~{g\/ F\) theor in the preasymptotic region used by some of the groups 
is not reasonable. 

In QCD, one can split g\ and Fi into leading (LT) and higher twist (HT) pieces 

g\ = 1?i)lt + (>i)ht, ^i = (^i)lt + (^i)ht- (1) 
Then, approximately 

gi _ (gl)LT ^ | (gl)HT (^l)HT j ^ 

^1 (F\)lt (gl)LT (F\)lt 

Note that LT pieces of g\ and F\ are expressed by the polarized and unpolarized PDFs, 
respectively. 

There are mainly two methods to fit the data - taking or NOT taking into account 
the HT corrections to g\. According to the first yj], the data on g\/F\ have been fitted 
taking into account the contribution of the first term h(x)/Q 2 in (^i)ht and using the 
experimental data for the unpolarized structure function F\ 
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According to the second approach [2] only the LT expression in (OQ) have been used in 
the fit to the g\/F\ data 
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It is obvious that the two methods are equivalent in a pure DIS region where HT 
can been ignored. To be equivalent in the preasymptotic region requires a cancellation 
between the ratios (gi )ht/(^i)lt and (F^ht/^Olt in ©• Fig. 1 (left), based on our 
recent results on g\ y|] and the unpolarized results of 04|], demonstrates the validity of 
this for x > 0.15, but clearly indicates that ignoring HT terms in the ratio gi/F\ below 
x — 0.15, as was done by groups using the second method or some its modifications 
[HI, @1, is incorrect. Note that in this case the HT effects are absorbed in the extracted 
PDFs which differ from those determined in the presence of HT (see Q3U for more 
details). 



POLARIZED PARTON DENSITIES AND HIGHER TWIST 

Let us start our discussion with the impact of the COMPASS [@] and CLAS [Q] inclusive 
DIS data on the longitudinal polarized PDFs and higher twist effects. The precise 
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FIGURE 1. Comparison of HT terms in g\ and F\ (left). HT effects for protons and neutrons (right). 




FIGURE 2. Impact of CLAS and COMPASS data on the uncertainties for NLO(MS) polarized PDFs. 



CLAS data on g*/Ff and g d jF( are data at low Q 2 (l -4 GeV 2 ) while COMPASS has 
presented data on Af mainly at large Q 2 and very small x (0.004 < x < 0.02), the only 
precise data at such small x. Compared to the HT values obtained in the LSS'05 analysis 
Q8Q , the uncertainties in the HT values at each x in the CLAS region (x ~ 0.1 — 0.6) 
are significantly reduced by the CLAS data, as seen in Fig. 1 (right). As expected, the 
central values of the polarized PDFs are practically not affected by the CLAS data. This 
is a consequence of the fact that at low Q 2 the deviation from logarithmic in Q 2 pQCD 
behaviour of gi is accounted for by the higher twist term in g\. The accuracy of the 
determination of polarized PDFs, however, is essentially improved. This improvement 
is illustrated in Fig. 2 (red curves). It is a consequence of the much better determination 




FIGURE 3. The three different LSS'06 solutions for AG(x) (left). Comparison with directly measured 
AG(x)/G(x) (right). 



of higher twist contribution in g\, as discussed above. The impact of COMPASS data 
on the PDFs uncertainties is also shown in Fig. 2 (the blue curves). As seen, they help 
to improve in addition the accuracy of the determination of the gluon and strange sea 
quark polarized densities at x < 0.2 and x < 0.1, respectively. Another important effect 
of COMPASS data (see |3|] and slides of Stamenov's talk at this Symposium Q9Q) is that 
As(x) and respectively, AL(x), are less negative at x smaller than 0.1. 

The present inclusive DIS data cannot distinguish between positive, negative and 
sign-changing AG(x) (see Fig. 3 (left)). In all the cases the magnitude of AG is small: 
|AG| < 0.4 and the corresponding polarized quark densities are very close to each 
other. COMPASS analysis finds also a negative AG(x), but has some peculiarities which 
suggest it is not physical. Also, we have found that if in the input gluon parametrization 
a term (1 + yx) is added, where 7 is a free parameter to be determined from the fit to 
the data, a negative AG solution is impossible. So, although the negative AG cannot be 
formally ruled out by the data, it seems to be not reasonable. 

A primary goal of the RHIC spin program is to determine the gluon polarization. 
There are two sources of information: p + p — > tz q (tz + I~) +X, p + p — > jet +X. We 
will mention here the combined QCD analysis of the inclusive DIS and RHIC 7T°- 
production data ft 1 Oil which has been recently performed by AAC 111 ill . Two solutions 
for AG(x), positive and node, have been found with values for the first moments at 
Q 2 = 1 GeV 2 : AG = 0.40 +/- 0.28 and -0.12 +/- 1.78, respectively. Note that for the node 
AG(x) solution the positivity condition AG(x) < G(x) is broken for x > 0.4, which makes 
it questionable. For more details about the possibilities of the RHIC data to constrain the 
gluon polarization AG see the talks presented by STAR and PHENIX Collaborations at 
this Symposium [9]. The main conclusion is that the magnitude of AG is smaller than 
0.4, but the present RHIC data cannot determine the form of the gluon density AG(x) . 

In Fig. 3(right) we compare the polarized gluon densities obtained by different QCD 
analyses (divided by the MRST'02 version of the unpolarized gluon density G(x)) with 
directly measured values of AG(x) / G(x) at [X 2 = 3 GeV 2 . One can see from Fig. 3(right) 
that the most precise values of AG/ G, the COMPASS ones |12[]. are well consistent with 




FIGURE 4. g p x at different Q 2 calculated using the LSS'06 PDFs corresponding to AG > and AG < 0. 

most of the polarized gluon densities determined in the recent QCD analysis, and one 
conclude that the present high p t measurements cannot distinguish between the different 
scenarios for AG(x) . 

It seems clear that present day data cannot distinguish between the three scenarios 
for AG(x). A clean distinction would be possible in the Electron-Ion Collider lfl3ll 
where very precise measurements of g^(x, Q 2 ) at very small x (x > 0.00075) could be 
performed. Fig. 4 shows gf (x, Q 2 ), calculated using the LSS'06 PDFs corresponding to 
positive and negative AG(x). There is a dramatic difference at small x. 

Very recently the DSSV group has presented results on polarized PDFs [14] obtained 
from the first global analysis of the polarized DIS, SIDIS and RHIC polarized pp 
scattering data. Due to the SIDIS data a flavor decomposition of the polarized sea is 
achieved. This analysis yields changing in sign As(x, Q 2 ): positive for x > 0.03 and 
negative for small x (see Fig. 5). Its first moment is negative (fixed in practice by the 
SU(3) symmetric value of a%) and almost identical with that obtained in the inclusive 
DIS analysis. It was shown in the talk by R. Windmolders at this Symposium [9J] that 
the determination of As{x) from SIDIS strongly depends on the fragmentation functions 
(FFs) and that for such an unexpected behavior of As(x) the new FFs |l5|l are crucial. 
So, the model independent extraction of FFs is very important. The NLO(MS) PDFs 
determined from analyses of different sets of data: LSS'06 (inclusive DIS data), AAC'08 
(iclusive DIS and RHIC ^-production data), DSSV (DIS, SIDIS and RHIC data) are 
compared in Fig. 5. Note that in the LSS analysis the HT corrections have been taken 
into account while AAC'08 and DSSV have used the second method described above 
(without account for HT terms). Although the first moments are almost identical, the 
quark densities themselves are different, especially As(x). Note that all analyses of 
the inclusive DIS data yield As(x) negative. So, the DSSV result on As(x) is a big 
challenge for our understanding of spin properties of the nucleon. We would like to 
emphasize also that the very accurate DIS data in the preasymptotic region require 
a more precise confrontation of QCD to the data, as was already mentioned in the 
Introduction. Otherwise, the results on PDFs will differ. A comment on an incorrect 
confrontation of QCD to the inclusive DIS data in the DSSV analysis is presented in 
Stamenov's talk [9] at this Meeting. 
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FIGURE 5. Comparison between LSS'06, AAC'08 and DSSV NLO PDFs in (MS) scheme. 



Let us finally discuss the present status of the proton spin sum rule. Using the values 
for AE(2 2 ) and AG(Q 2 ) at Q 2 = 4 GeV 2 obtained in LSS'06 analysis [H one can find 
for the spin of the proton (the numbers in brackets correspond to node AG): 

S z = l - = ^AE(£ 2 ) + AG(£ 2 ) +L Z (Q 2 ) = 0.55(0.15) ±0.25(0.49) +L Z (Q 2 ). (5) 

Although the central values of parton contribution are very different in the two cases, in 
view of the big uncertainty in © coming mainly from the gluons, one cannot make a 
definite conclusion about the quark-gluon contribution in the spin of the nucleon. 
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